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ABSTRACT 



Context. According to the Big Bang theory and as a consequence of adiabatic expansion of the Universe, the temperature 
of the cosmic microwave background (CMB) increases linearly with redshift. This relation is, however, poorly explored, 
and detection of any deviation would directly lead to (astro-)physics beyond the standard model. 

Aims. We aim at measuring the temperature of the CMB with an accuracy of a few percent at 2=0.89 toward the 
molecular absorber in the galaxy lensing the quasar PKS 1830—211. 

Methods. We adopt a Monte-Carlo Markov Chain approach, coupled with predictions from the non-LTE radiative 
transfer code RADEX, to solve the excitation of a set of various molecular species directly from their spectra. 
Results. We determine Tcmb=5.08±0.10 K at 68% confidence level. Our measurement is consistent with the value 
rcMB=5.14 K predicted by the standard cosmological model with adiabatic expansion of the Universe. This is the most 
precise determination of Tcmb at z>0 to date. 



Key words, cosmology: observations - cosmic background radiation 
lines - quasars: individual: PKS 1830—211 



cosmological parameters - quasars: absorption 



1. Introduction 

1.1. The CMB temperature as a function of redshift 

The cosmic microwave background (CMB) is one of the pil- 
lars of the Big Bang theory. This radiation field is a relic of 
the epoch when matter and photons decoupled at z^^llOO 
and fills up the entire sky with a photon density of ^400 
cm~^. Since decoupling, the CMB photons have cooled 
down with the expansion of the Universe. Observations with 
the COBE satellite have demonstrated that the CMB cor- 
responds to a nearly perfect black body, characterized by a 
temperature Tg at z=0, which is m easured with very high 
accuracy, To=2. 72548 ± 0.00057 K ()Fixsenll2009t ). 

Measuring the CMB temperature at high redshift has 
considerable cosmological interests, in 1) demonstrating 
that the CMB radiation is universal (Equivalence principle) 
and 2) tracing the evolution of its temperature with red- 
shift, Tcmb{z)- Adiabatic expansion predicts that the CMB 
temperature evolution is proportional to (l+z). Alternative 
cosmologies, such as decay i ng dark energy mode ls {e.g. 
iLimalllQQdiLima et al.ir2000HJetzer fc T ortora"201l') where 
dark energy can interact with matter via creation of pho- 
tons, affecting the CMB spectrum, predict a deviation from 
this simple law with T(z)=Tqx(1+z)^^~°'\ Precise mea- 



surements of Tcmb at high redshift can therefore constrain 
those theories. 

Two methods can be used to probe the CMB temper- 
ature at z>0. The first one is based on multi-frequency 
Sunyaev-Zel dovich (S-Z) observat ions toward galaxy clus- 
ters (see e.g. iHorellou et aLll2005[ ). Inverse Compton scat- 
tering by free electrons in a cluster of galaxies induces a 
distortion of the CMB spectrum, which can be used to 
con strain the CMB t e mpera ture at the location of the clus- 
ter. iBattistelli et al.l ()2002[ ) reported the first determina- 
tion of Tcmb{z) toward the Coma (z=0.02) and Abell 2163 
(z=:0.2) clusters. More recently, this method was applied 
to a sample of 13 clusters between z =0.023 and z=0.55 
by iLuzzi et all ([200-9), who found 0=0.024155:°^^, thus a 
TcMB~z law in agreement with the standard model. One 
limitation of the S-Z m ethod is the scar c ity of clusters 
at high-redshift iz>l). Ide Martino et al.l (|2012l) discuss 
prospects in the near future based on Planck data. 

The second method relies on spectroscopic studies 
of lines in absorption against quasars and their excita- 
tion analysis. Most of the measurements at high red- 
shi ft have used UV spectroscopy of ato r nic species (Ci : 
Meyer et all Il986t ISongaila et al.l Il994b|: lGe"et all 11997 

ID I 



iRoth fc Baueilll999t Cii: ISongaila et al.[ll994al: iLu et all 
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119961: iMolaro et al.ll2002l: Ci & C ii: ISrianand et al.ll2n0(l . 

Strictly speaking, these measurements are all upper limits 
on Tcmb, as the contributions from other local sources of 
excitation (collisions, local radiation field) are yet largely 
uncertain and poo rly constrained, and have t o be accounted 
for. More recently. iNoterdaeme et all ()2011[ ) observed elec- 
tronic transitions of CO absorption toward several quasars 
at redshift between 1.7 and 2.7. Together with the previ- 
ous measurements of Tcmb(-z) in the literature, they could 
obtain the tightest constraint to date on the temperature 
evolution of the CMB, a=-0.007±0.027, consistent with 
adiabatic expansion of the Universe. Nevertheless, all the 
studies so far have relied on a limited number of transi- 
tions/species, introducing a possible bias in the determi- 
nation of excitation conditions. Also, the UV spectroscopy 
technique {e.g. of Ci and Cii) is not suitable below Z'^1, 
due to the atmospheric cut-off around 300 nm. 

The use of radio- mm molecular absorbers to determine 
Tcmb at z>0 is another particularly attractive method, 
due to the tighter constraints on physical conditions and 
local excitation when a variety of molecular species is de- 
tected. In this paper, we follow this approach, with the aim 
to determine the CMB temperature from multi-transition 
multi-species excitation analysis in the z=0.89 radio-mm 
molecular-rich absorber toward the quasar PKS 1830—211. 

1.2. The z—0.89 molecular-rich absorber toward 
PKS 1830-211 

PKS 1830— 211 is a radio- loud blazar at a redshift z=2.5 
(|Lidman et al.llT999l ). gravitatio nally lensed by a z=0 . 88582 
almost face-on s piral galaxy (jWiklind fc Combei 119961 : 
|W mn et all I2002D . The presence of an additional galaxy 
at z= 0.19 on the line of sight was inferred from H i absorp- 
tion ()Lovell et al.lll996D . The lensed image of the quasar 
appears in the radio continuum as a textbook example of 
a complete Einstein ring of 1" in diameter, plus two bright 
and c ompact com ponents, located to the NE and SW of the 
ring (jJauncev et al. 1991). The Einstein ring has a steep 
spectral index and grows fainter with frequency, until in- 
visible in the millimeter. The NE and SW images of the 
quasar are marginally resolved in 43 GHz VLBI obse rva- 
tions with 0.5 mas beam resolution bv lJin et al.l (|2003[ ). At 
the distance of the z=0.89 galaxy, their angular size cor- 
responds to a scale of ~1 pc at mm wavelengths. Most in- 
terestingly, the line of sight toward both images intercepts 
molecular clouds located in the spiral arms of the inter- 
vening z^O.89 galaxy, which results in remarkable molec- 
ular absorptions ( W ikli nd &: Combes 1998; Menten et al.l 
[1991 iMuller et al.M 2006D. In fact. PKS 1830-211 is the 
strongest of on ly fou r redshifted mm-molecular absorbers 
known to date ((Com bes 2008), the one at th e highest red- 
shift, and the most molecule-rich (jMuller et al.l 1201 1[ ). As 
such, PKS 1830—211 is a target of choice for using molecules 
as cosmological probes, for example of the CMB tempera- 
ture, of the constancy of fundamental constants , or of the 
chemical evo lution of the Universe (see e.g. Henk el et al.l 
r2009t IMulleill20Tll) . 

The activity of the blazar leads to time variations of its 
flux, with a time delay between the NE and SW images 
of ab out 25 days (iLov ell et al] 119981: I Wiklind fc Comljesl 
Il999i) . as well as changes in the morphology and barycen- 
ter of the continuum emission, seen in VLBI observat ions 
(|Garrett et al.lll997t iJin et al.l [20031 ). iNair et al.l (12005D in- 



terpret the changes in the apparent distance between the 
NE and SW images by the recurrent ejection of plasmons 
along a helical jet in the blazar. Consequently, and of great 
importance for our study, the continuum emission of the 
two images is drifting behind the abso r bing gas, at the 
rhythm of plasmon ejections. iJin et al] (|2003D measured 
drifts up to 200 /ias (1.6 pc in the lens plane at z=0.89) over 
eight months. The changes of sightline result in a time vari- 
ation of the molecular absorption profiles, as revealed by 
the long-term monitori ng of the HCO"''/IICN 1-0 lines by 
iMuller &: GuelinI (|2008[) . This time variation could hamper 
studies which rely on the comparison of different molecular 
transitions, such as constraints on the constancy of fun- 
damental constants, multi-transition excitation analysis, as 
well as abundance and ratio measurements. 

2. Observations and data reduction 

Our main goal was to observe several rotational transi- 
tions of various molecular species in order to derive their 
excitation and constrain the CMB temperature toward 
PKS 1830—211. In particular, to take advantage of the two 
SW/NE lines of sight through the intervening galaxy, which 
are independent (on either side of the bulge, at galactocen- 
tric distances of ~2 kpc and ~4 kpc, respectively) and likely 
to have different excitation conditions, giving a robust mea- 
surement of Tcmb at z=0.89. 

Toward the NE li ne of sight, only a handful of species 
have been detected (jMuller et al.l 1201 it) : we have chosen 
to focus on the HCO+, HCN, and HNC J=l-0, 2-1, and 
3-2 transitions, which are all optically thin. These same 
transitions arc saturated or optically thick toward the SW 
line of sight, but their isotopic variants, observable in the 
same tunings, can be used instead. Several other species, 
such as HOC+ and SiO, also have rotational transitions in 
the same observable bands. 

In order to cover the different rotational transitions 
of our targeted molecular species, we used the Australia 
Telescope Compact Array (ATCA) at 7 mm and 3 mm 
and the Plateau de Bure interferometer (PdBI) at 3 mm 
and 2 mm. Details on the respective observations and spe- 
cific data calibration are given hereafter in ^ 12.11 and ^ \2.2[ 
and a summary of the observations is given in Table [TJ 
Unfortunately, because of scheduling and weather con- 
straints the ATCA and PdBI observations could not be 
carried out simultaneously, which raises the concern of time 
variations of the source between the observations. This 
problem is discussed in ii l2.3l 

The full final spectra are shown in Fig lA.ll and IA.2[ A 
crude comparison of the new data and ATCA spectra ob- 
tained in 2009/2010 shows that significant variations have 
occurred fFig jA.2|) : while the intensities of the SW absorp- 
tion lines (ti=0 kms^^) have increased by a factor of ^2 
for most species, the absorption depths of the HCO"*" and 
HCN hncs toward the NE image (z)=-147 kms~^ El) have 
significantly decreased. The resulting low signal-to-noise ra- 
tio on these lines and the continuum baseline uncertainty 
due to the wings of the corresponding saturated lines from 
the SW component (see Fig |A.2[) prevent us from obtaining 
a robust fit, which hampers our initial goal to derive Tcmb 
in several independent lines of sight. 

^ All velocities throughout this work refer to the heliocentric 
reference frame, adopting a redshift z=0. 88582 
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Table 1. Summary of the observations. 



Date 


Telescope 


Band 


Frequency 


SV " 


a' 




R 


Notes 








(GHz) 


(kms-^) 


(10-^) 








2011 Jul. 24 


PdBI 


3 mm 


93.0 - 96.6 


-/2.0 " 


-/2.3 <= 


0.416 ± 0.006 


1.40 ± 0.03 


Total flux=2.6 Jy 


2011 Jul. 27 


ATCA 


7 mm 


44.8 - 48.9 


6.3 


2.0 


0.402 ± 0.003 


1.49 ± 0.01 




2011 Jul. 27 


ATCA 


3 mm 


90.8 - 96.4 


3.2 


2.8 


0.397 ± 0.004 


1.52 ± 0.02 




2011 Aug. 17 


PdBI 


2 mm 


140.2 - 143.8 


4.1/1.3 


4.5/6.8 " 


0.437 ± 0.010 


1.29 ± 0.05 


Total flux=2.1 Jy 



Notes. " Velocity resolution for the HCO''' lines. '' Noise level normalized to the total continuum (NE+SW images). Fraction of 
the total continuum intensity corresponding to the SW image, measured on the saturated part of the spectra between —7 kms^^ 
and +7 kms^^. '* Magnification factor R=/i,g(NE)//(,g(SW), assuming no flux comes from other components. ^ Corresponding to 
the WIDEX and narrow-band correlator output, respectively. 



Other velocity components (e.g. at vel ocity 
-224 kms^i and +170 kms-\ see lMuller et al.ll201lD . are 
also too weak for a meaningful analysis. 

2.1. Australia Telescope Compact Array (ATCA) 

Observations were carried out with the ATCA on 2011 
July 27, with four antennas in a compact configuration. 
The longest baseline was ~150 m, providing an angular 
resolution >9" at 7 mm and >4" at 3 mm. The contin- 
uum emission (1" Einstein ring and NE/SW images) of 
PKS 1830—211 was therefore not resolved. 

The C ABB correlator was set up to provide a bandwidth 
of 2 GHz with a spectral channel spacing of 1 MHz in each 
sideband. We used two partly overlapping tunings at 7 mm 
and three at 3 mm to cover the frequency ranges given in 
Table[T] The on-source integration time was 0.5 h/tuning at 
7 mm and 1.5 h/tuning at 3 mm. The total allocated time, 
including calibration and overheads, was 8 h. The receivers 
were tuned in dual polarisation mode. Both polarisations 
were averaged after the data calibration. The bandpass cal- 
ibration was performed on the bright quasar 1921—293 for 
each tuning. 

The data reduction was done using MIRIAD 
(jSault et al.l 119951 ). Source visibilities were self-calibrated 
and the (unresolved) continuum level normalized to unity. 
Spectra were then extracted from the calibrat ed visibilities 
following the same method as described by iMuUer et al.l 
((2011). A fixed Doppler track correction of +13 kms~^ was 
applied to all data. The drift in velocity being less than 
0.1 kms~^/hour, the effect introduced by adopting a fixed 
correction is negligible considering the spectral resolution. 

In complement to this dataset, we use previous ATCA 
7 mm o bservations from 2 009 September 1 and 2 pub- 
lished in lMuller et al.l (|201l[ ). specifically for transitions of 
C-C3H2, CH3CN, SO, and HC3N (see §EXl|). These ob- 
servations were obtained and reduced following the same 
procedure as described above. 

2.2. Plateau de Bure interferometer (PdBI) 

Observations were obtained with the PdBI on 2011 July 
24 at 3 mm, and on 2011 August 17 at 2 mm, with five 
antennas in a compact array. All baselines were shorter than 
100 m, so that the continuum emission of PKS 1830— 211 
was not resolved, as for ATCA observations. 

^ http:/ /www. atnf.csiro.au/computing/software/miriad/ 



The wide-band correlator WIDEX was used, and pro- 
vided a 3.6 GHz bandwidth with a spectral resolution of 
1.95 MHz. In addition, the narrow-band correlator was con- 
figured to select two 160 MHz spectral bands with a bet- 
ter spectral resolution of 0.625 MHz, with dual polarisa- 
tion. At 3 mm, these narrow bands were set to cover the 
kms^^ and —147 kms~^ absorption lines of both HCO+ 
and HCN 2-1. At 2 mm, the narrow-band IF range allowed 
us to cover only the kms~^ and —147 kms~^ components 
of the HCO"'" 3-2 line and the —147 kms^^ component of 
the HCN 3-2 hne (see the spectra in FigH]). 

The bandpass was calibrated on the bright quasars 
1749-1-096 (at 3 mm) and 3C454.3 (at 2 mm). We fol- 
lowed the standard calibration procedure in CLIC 0, using 
the "Self-Calibration on point source" routine. The spec- 
tra were then exported to CLASS ^, where we divided the 
spectral baseline with polynomials of degree two for the 
narrow-band spectra, and degree 12 or 24 depending on 
the case for the full 3.6 GHz WIDEX bands, to correct for 
residual ripples in the spectral baseline. The target absorp- 
tion lines are too narrow to be seriously affected by these 
polynomial orders. 

Features corresponding to O3 atmospheric lines were 
identified in the final spectra. In particular, we flagged a 
region of 200 MHz around the 142.175 GHz O3 line in the 
2 mm data. 

2.3. Continuum illumination 

The absorption intensities (lu), measured from the total 
continuum level, are converted to optical depths {t^) ac- 
cording to: 

where Ibg and fc are the fraction of the total continuum 
intensity (normalized to unity) and the continuum source 
covering factor, corresponding to each individual image 
of the quasar, respectively. We adopt a value /c=l (see 
iWiklind & Combes 19J8|; iMuller fc Gu"Ih] [2008) . 

At radio wavelengths, the continuum emission of the 
background quasar is dominated by two bright and compact 
images to the NE and to the SW of the 1" Einstein ring (see 
iJ1.2p . These NE and S W images are not resolved in our ob- 
servations, so that we cannot measure directly their respec- 
tive fluxes /bg(NE) and /bg(SW). However, the HCO+ and 

^ http://www.iram.fr/lRAMFR/GlLDAS/ 
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HCN lines appear to be flat-bottomed near v~0 kms^^, 
and both at the same intensity for a given epoch (FigllJ, 
which strongly suggests that they are saturated. The de- 
tection of the ^^C variants of these species (and even 
of other less abundant isotopologues) clearly implies that 
they are optically thick, with opacities higher than t=10. 
Because we know that the absorption near v=0 kms~^ oc- 
curs in front of the SW image (Wiklind fc Combes! [19981 
iMuUer et alJ l200l . we can measure Ibg{SW) from the 
saturation level of the HCO"'" and HCN lines between 
v=—7 kms""'^ and v=+7 kms^^, under the assumption 
fc=l- The results are given in Table [T] The comparison of 
the HCO+ J=3-2 spectra obtained with the WIDEX cor- 
relator (velocity resolution SV of ~4 kms^^) and with the 
narrow-band correlator {6V^1 kms~^) shows that the lim- 
ited resolution in velocity does not affect significantly the 
determination of the saturation level. 

The measured values of /bg(SW) differ between the 
three epochs, with the biggest change on the last measure- 
ment with the PdBI, obtained 20 days after the ATCA 
observations. As the saturation level does not depend on 
the absolute flux density calibration (the intensities are 
normalized to the total continuum level), we believe that 
these changes are significant, and not of instrumental ori- 
gin. They could be caused either by a change of the NE/SW 
ratio due to the tim e delay of '-^25 days between the two im- 
ages of the quasar (iLovell et al.l [19981: iWiklind fc Combei 
[1999*), or to a change of the source covering factor if the 
sightline moved between the observations. Changes of the 
continuum morphology, related to the jet activity of the 
blazar, are clea rly seen within a time span of two weeks 
(jJin et al.ll200l see also lGarrett et al.lll997D . but monitor- 
ing on a shorter timescale is still lacking, to the best of our 
knowledge. 

The changes seen between July and August 2011 are 
unfortunate for our purpose, because the lines of sight 
through the intervening galaxy - hence the column den- 
sity and physical conditions of the absorbing gas - might 
have changed. This emphasizes the need to conduct the ob- 
servations within a short time interval when the goal is to 
combine different observations toward PKS 1830—211. In 
the rest of the paper, we will treat only the ATCA observa- 
tions of the lines corresponding to the SW component, and 
assume that no significant intra-day variations occur. 

In order to take into account the uncertainty on the 
continuum illumination A/hg, we have added quadratically 
to the fit uncertainties of the integrated opacities a term 
calculated as: 



rdv 



rdv 



X 4^ X Fil/h,), 



ho 



(2) 



where the function F{x) = x/{{x — l)ln(l — x)), as ob- 
tained from the differentiation of EqlT]with respect to Ibg. 
This correction is simply equivalent to Al — for opti- 
cally thin lines. 

The magnification ratios R=NE/SW, which we take 
as the fraction (1— /(,g(SW))//f,g(SW), are listed in 
Ta ble [J They are within the range of values compiled 
bv iMufler fc GuelinI (l2008h . using the same measurement 



method of the saturation level of the HCO+ SW absorp- 
tion. 



3. Analysis 

In this section we describe the analysis of the spectral data, 
starting from a simple approach to derive rotation temper- 
atures. We then investigate the effects of local physical con- 
ditions on the line excitation. Finally, we provide a global 
solution of the problem using a Monte-Carlo Markov Chain 
method, from which we derive our final measurement of 
Tcmb at z=0.89. 



3.1. Rotation temperatures 

First, we focus on deriving the rotation temperature of the 
observed species assuming that the excitation is dominated 
by radiative coupling with the CMB (low-density case). 
Under this condition, the population of the energy levels 
follows a Boltzmann distribution, for which the rotation 
temperature equals Tcmb- The rotation temperature Trot 
between two energy levels (upper u and lower I) is defined 
by the Boltzmann equation: 



— = — exp -TTjr- 

ni gi \ fcTi-ot 



(3) 



where and ni are the populations, gu and gi the degen- 
eracies, of the upper and lower level, respectively, connected 
by a transition at frequency Vui ■ The column density A/'lteB 
can then be derived as: 



A^LTE — 



3h Q(^™')^^p(^Jt;) 



1 — exp 



V ksl rot / 



rdv, 



(4) 



where Q(Trot) is the partition function, Ei the energy of 
the lower level with respect to ground state, the dipole 
moment, v is the frequency, Sui the line strength and J rdv 
the integrated opacity of the line. 

The kinetic temperature of the gas along the SW line 
of sight is of the order of ~80 K, as indicated by obser- 
vations of the symmetric top molecule s NH3 and CH3CN 
(|Henkel et al.l 120081 : [Miller et al.l 1201 ill . Bec ause the den- 
sity is moderate (of the order of 10^ cm~^, iHenkel et al.l 
2009), the gas is sub-thermally excited and the excitation 
mostly coupled radiatively with the CMB. Hence, the rota- 
tion temperatures are expected to be close to (only slightly 
higher than, see N3.2p Tcmb- The ro tation temperatures 
of sev er al species have b een d er ived by Combes Sz ^yiklindl 



(I1999D. iMenten et al.l ([19990 . iHenkel et al.l (|2009f) . and 



iMuUer et al.l (|201lD in the SW line of sight, and are consis 



tent with this picture. 

We note that the Ray leigh- Jeans approximation does 
not hold for v>100 GHz transitions when the temperature 
is low (~5 K), and t he rotation diagram method (as done 
in lMuUer et al.|[2011l for lower frequency data) is not appli- 
cable to derive the rotation temperature here. Instead, we 
use EqHl and plot the A^LTE^Trot curves for the different 
observed transitions (Figl^)). These curves have slightly dif- 
ferent sensitivities to the temperature. Their intersection 
gives the rotation temperature and column density. The 
results are given in Table [51 As expected, the rotation tem- 
peratures are all close to ^5 K. 



* We refer to this estimate as the LTE column density, al- 
though strictly speaking, we simply assume that the excitation 
is dominated by the CMB photons. 
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Table 2. Results for rotation temperature measurements 
toward the SW absorption. 



Species 


Dipole 


Date 


Trot 




moment 


or the 


(K) 




(Debye) 


observations 


C2H 


0.77 


2011 


5.3 ± 0.1 


bO 


1.54 


2009 


5.4 ± 1.4 


HNCO 


1.58 


2011 


9.8 ± 1.5 






901 1 


0.1 m u.-^ 


H"CN 


2.99 


2011 


5.1 ± 0.2 




2.99 


2011 


4.1 ± 0.4 


HNC 


3.05 


2011 


4.6 ± 0.2 t 


HN"C 


3.05 


2011 


4.8 ± 0.3 


SiO 


3.10 


2011 


6.0 ± 0.2 


C-C3H2-0 


3.43 


2009 


5.6 ± 0.4 


C-C3H2-P 


3.43 


2009 


5.4 ± 1.0 


HC3N 


3.73 


2009 


6.3 ± 1.3 


H"CO+ 


3.90 


2011 


5.3 ± 0.1 



Notes. ' Might be affected by opacity effects. Uncertainties are 
given at la. 



3.2. Effects of local physical conditions 

In order to test the impact of the local physical condi- 
tions on the line excitation, we have run models using 
the non-LTE molecular radiative transfer code RADEX 
(|van der Tak et all l2007fl . varying parameters such as the 
kinetic temperature (Tkin), the H2 density (^Ha), and the 
local radiation field. 

RADEX uses the escape probability method for a ho- 
mogeneous medium and is based on given radiative and col- 
lisional transition rates, which are only available for a few 
species in a limited range of temperatures and column den- 
sities. In this work, we use the standard files available from 
the Leiden Atomic and Molecular Database (LAMDA , 
ISchoier et al.l [20051 ) . except for IIC3N, for which rates ex- 
tended to 100 K we re kindly provided to us by A. Faure 
(jWernh et al.ll2007allbl . and Faure & Wiesenfeld m prepara- 
tion) . We consider that H2 molecules are the main collision 
partners. 

3.2.1. H2 density and kinetic temperature 

We have run RADEX simulations varying riHs between 
500 and 10^ cm~^ for two representative values of the 
kinetic temperature Tkin=50 K and 100 K. These ranges 
of values enc ompass previous e stimates of the density 
(~2000 cm-3,'Henkel et al."2009^, and of the kinetic tem- 
perature (~80 K, Hcnkcl ct al. 2008) toward the SW com- 
ponent of PKS 1830—211. Such a value of the kinetic tem- 
perature is typical of T\\-n observed in Galactic diffuse 
clouds {e.g. ISnow fc McCalll [2006). At this stage and for 
sake of simplicity, we have set the column densities of the 
various species to the values that can be derived directly 
from Fig 12] from the intersection of the A^lte vs T^ot curves. 

The results of our RADEX simulations can be seen in 
FiglSl The curves plotted on this figure indicate the evo- 
lution of the rotation temperature for a given transition 
as a function of the H2 density, for HCO+, HCN, HNC, 

^ http:/ /home. strw.leidenuniv.nl/~moldata/ 



some of their isotopologues, HC3N, SO, HNCO, SiO, and 
C-C3H2 (ortho and para forms). The rotation temperature 
measured from the observation of two transitions is then 
the average of the corresponding curves for these two tran- 
sitions. The rotation temperature increases with H2 density, 
as can be expected when collisional excitation progressively 
adds to the radiative excitation from the CMB. This is the 
case of all transitions, except for the c-C3H2(para) 331-322, 
which shows a decrease of Tiot with njjj • 

Among the various species investigated here, the most 
sensitive to collisional excitation are HC3N and HNCO, 
for which AT=Trot— ?cmb reaches values >1 K for 
nH2>2000 cm~'^. On the other hand, the transitions of 
HCN, HNC and C-C3H2 (ortho) (and their isotopic variants 
when applicable) have AT of less than a few tenths of a 
Kelvin for riHa'^SOGO cm~^, and can therefore be consid- 
ered as good Tcmb indicators. 

3.2.2. Local radiation field 

We now add the excitation from a local radiation field. We 
consider first the mean background radiation near the Sun's 
location in the Galaxy, which consists mostly in dust ther- 
mal emission (we adopt a single-tem perature component 
with Td„st=23 K. IWright et al.lll99lD . Non-thermal radia- 
tion is negligible at mm/submm wavelengths. We did not 
find measurable differences as compared to the previous 
case with CMB excitation only. Gradually increasing the 
dust temperature, we start to see some effects on the rota- 
tion temperature of HNCO for Tdust^lOO K, while all other 
molecules do not show any changes. We thus conclude that 
the effect of the local radiation on the excitation of the 
observed lines is negligible, unless very peculiar conditions 
prevail. 

3.3. A global solution for the excitation conditions with 
Monte-Carlo Markov Chains 

In the following, we adopt a global method for solving si- 
multaneously the physical conditions and Tcmb- We first 
describe our working assumptions and the RADEX mod- 
eling, and then the Monte-Carlo Markov Chain approach 
used to adjust the model to the spectra. 

3.3.1. RADEX model 

For a given molecular species, we use RADEX to com- 
pute the velocity integrated intensity of each line present in 
our spectra, given the background radiation field (Tcmb), 
the kinetic temperature (Tkin), the number density of H2 
molecules (riHa), the column density of each molecule (iVcoi) 
and the width of the molecular line (AV), which is assumed 
to be Gaussian. 

The physical conditions are possibly changing with time 
(whenever the barycenter of the continuum emission drifts). 
We therefore consider independent physical condition sets 
for both the 2009 and 2011 data. For simplicity, we as- 
sume that the absorbing gas is homogeneous and can be 
characterized, for each epochs, by a single specific kinetic 
temperature Tkin and density, common for all molecu- 
lar species. For the SW line of sight, we thus have the set of 
parameters Tkin/»^H2(2009) and Tkin/»^H2 (2011). The back- 
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ground radiation field is assumed to be constant and equal 
to the CMB. 

Given the velocity centroid (Vq) and linewidth (AV), 
specific for each molecular species, we are then able to con- 
struct a synthetic spectrum, using Eq. (TJ covering the ob- 
served frequencies. For each species, we thus have three 
free parameters (-/Vcoi, AV, and Vq) and three tied param- 
eters (Tkin, and Tcmb) to fit. In the fit, we include 
all species with molecular data available in the LAMDA 
database, and with at least two, non saturated, lines in 
the observed spectra. We end up using: H^'^CO+, H^-^CN, 
HC^^N, C-C3H2 (both ortho and para forms), SiO, CH3CN, 
SO, HC3N, and HNCO. 

3.3.2. Adjusting the model 

Such a high number of variables is very difficult to constrain 
using classical techniques. Our approach here is to adopt a 
MCMC meth od, widely used to es timate cosmological pa- 
rameters (e.g. tounklev et al.l 12005' and references therein). 
We use a subclass of the MCMC methods, the Metropolis- 
Hastings algorithm to probe the probability distribution of 
our problem and study the degeneracies between the differ- 
ent parameters. 

For each parameter, we choose the proposed density, or 
jumping distribution, as Gaussian distributions, with vari- 
ances set initially to the result of the estimate made in 
Muller et al. (2011). As the collisional rates are only tab- 
ulated over a particular temperature and density range, 
we rejected all jumps going outside the allowed range. The 
likelihood is simply computed by p oc exp (— x^/2), with 



N , 



od) 



i/=l 



(5) 



where Sy is the observed spectrum with its uncertainties 
(7 J/ at frequency channel and Si^^mod the modeled spectra 
with its uncertainty (T^^mod in order to take into account 
the uncertainty on the continuum illumination A/f,g (see 
EqiD. 

After a few initial chains, the variances were set ac- 
cording to the results of the chains themselves, in order 
to speed up the convergence of the adjustment. Note that 
the final result does not depend on the proposed density 
distribution, as long as it is symmetric and the chains 
ha ve converged, ^ ye us ed the spectral analysis described 
in iDunklev et al.l ()2005() to ensure the convergence of our 
chains before doing the analysis. 

3.3.3. Result of the Monte-Carlo Markov Chain method 

One major advantage of the MCMC approach is that it al- 
lows us to take all uncertainties into account when deriving 
the posterior distribution of the parameters. Most of the 
parameters display a posterior Gaussian distribution with 
median value and 68 percentile error reported in Table |31 
Fig. m presents the results for a subset of the parameters, 
Tcmb and, for the two observing campaigns, Tkin and riHj- 
As expected, there is a slight degeneracy between Tcmb 
and njja , for both campaigns, where a higher background 
radiation field could be compensated by a lower H2 density. 
However, the number of lines used here allows us to remove 
partially this degeneracy, and to put firm constraints on 



Tcmb- A percentile analysis (68% confidence level) of the 
final chains, including all lines from all species for both 
campaigns, leads to a value of the CMB temperature of 



Tcmb{z = 0.89) = 5.0791 



+0.0993 
0.0994 



(6) 



4. Discussion 

4.1. Cosmological constraints 

Our measurement of Tcmb at z=0.89 is interesting because 
it is the most precise measurement at z>0 to date, reach- 
ing a remarkable precision of 2%; ii) the observations of 
various molecular species allow us to constrain the physical 
conditions of the absorbing gas, and limit the possible bias 
of using only one species; and in) it is derived at a redshift 
intermediate between S-Z measurements in galaxy clusters 
{z<Q.6 so far) and optical/UV quasar absorption systems 
{z>l.b due to atmospheric cut-off), and with an alternative 
method (molecular excitation). 

Including our value of the CMB temperature at 
z=0. 89 with other re s ults published in the literature 
(see iNoter daeme et all I2OIII and references therein for 
points plotted in FiglSJ, we obtain a tighter constraint 
on the Tcmb(^)=Tox(1-Hz)(i-") law, a=-h0.009±0.019, 
compared to the constra int a=— G.007±0.027 reported by 
INoterdaeme et all ()201lD . As far as we can tell, the CMB 
temperature evolution is consistent with the standard cos- 
mology model, with adiabatic expansion and a—Q. 
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Fig. 5. Measurements of the CMB temperature as a func- 
tion of redshift. Data points in green correspond to S- 
Z measurements toward galaxy clusters, in black to C i 
and C II absorption stud ies, in blue to CO absorption (see 
INoterdaeme et al] 1201 ll and references therein) , and the 
value derived toward the PKS 1830— 211 SW absorption 
is marked in red. The dotted line corresponds to the law 

TcMB=ToX(l+z). 



4.2. Physical conditions of the absorbing gas toward the SW 
component 

Besides Tcmb, our MCMC approach also allows us to 
constrain the physical conditions toward the SW compo- 
nent, independently for the 2009 and 2011 observations 
(Table 121). We find that the kinetic temperature is well 
constrained, mostly owing to the observations of different 
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Table 3. Results (median and 68 percentile error) of the Metropolis-Hasting estimates of the physical conditions toward 
the SW component of the 2009 and 2011 campaigns. 



Species 


Tcmb 




Tkin 


log TVcoi 


AV 


Vo 




[K] 


[cm-^] 


[K] 


[cm-2] 


[kms-i] 


[kms-i] 


-0.10 



+0.3 
-0.3 
+0.5 
-0.5 



(2009) - 800 ^400 81 

CH3CN - - - ' 13.200 tool? 20.0 tol -2.3 tot 

SO _ _ _ 13.458 19.8 -1.8 

c-C3H2(o) _ _ _ 13.711 tHle 24.1 Ig-.^ -3.1 1°:^ 

c-C3H2(p) - - - 13.223 tHlg 21.0 t°ol "l-S 

HC3N - - 13.037 22.3 t^;^ -5.6 

(2011) - 2200 tlZ 82T? 

CH3CN _ _ _ 13.530 +g 024 Q +1.1 _3 9 +0.5 

_ _ _ 13.264 ™ 16.2 Ig-^ -2.4 l^^i 

Hi^CN _ _ _ 13.385 17.0 1°-^ -2.8 l^'^ 

HC^^N _ _ _ 12.607 ™ 14.2 1^8 -3.2 l^'^ 

HNCO _ _ _ 13.662 t° °^9 ^4 5 +1.2 q +o.4 

SiO - 13.369 ™ 14.6 1^3 -3-0 Ig:! 

Notes. A^coi, AV, and Vo are the column densities, linewidths (FWHM), and velocity centroids, respectively. 



K-transitions of CH3CN. Our result of Tkin'^80 K is compa- 
rable for both epochs, and also similar to the value derived 
bv lHenkel et al.l ( 2008) from observations of NH3 inversion 
transitions in 2003-2004. Despite probable changes in the 
line of sight, it seems that the kinetic temperature of the 
absorbing gas remains unchanged. The situation is not so 
clear for the H2 volume density, whic h is not as we l l con - 
strained as Tki„. From 2001-2002 data. lHenkel et al.l (|2009D 
estimate nH2~2000 cm~'^, comparable to our 2011 result. 
By comparison, n-n^ for our 2009 data is roughly half that 
value, a result that is also found for the integrated opaci- 
ties (or column density, all other parameters being equal) 
between 2011 and 2009. It is possible that the 2011 line 
of sight has intercepted a denser, more compact, absorbing 
component than in 2009. In any case, based on the observed 
rotation temperatures and the value of Tkin, it is clear that 
the excitation of the absorbing gas is sub-thermal, and 
that the density should be moderate (<10'* cm""^). This 
is consistent with the absorbing gas being composed of dif- 
fuse/translucent clouds (see .Muller et al. 2011.). 

We note that the increase of absorption intensities in 
2011 allows us to detect several new species toward this 
galaxy (Table [Ql: NH2CHO (2 transitions), ^OSiO (J=2- 
1), HCS"*" (whose J=2-l transition shows up better out of 
the w ing of the C-C3H2 2i2-loi hue than in the lMuUer et al.l 
I2OIII data), and tentatively HOCO"*" (detected at just Aa 
at rest frequency 85.531512 GHz). 

4.3. Limits and Perspectives 

Our excitation analysis relies on the simplifying assumption 
of an homogeneous medium, characterized by single values 
of the H2 density and kinetic temperature common for all 
molecules, when it is an observational fact that Galactic 
clouds are highly structured and turbulent. iHenkel et al.l 
(|200l) interpret NH3 observations to show two components 
of absorbing gas, one at approximately 80 K occupying 
80-90% of the column, and a second at kinetic tempera- 
ture greater than 100 K occupying 10-20% of the column. 
The warmer component would be even more diffuse and 
more subthermally excited, if in pressure equilibrium with 



the denser and cooler phase. The typical scale length cor- 
rcspond ing to a H2 column density of 2x10^^ cm~^ (see 
§4.2 bv lMtdier et aLllMll for the SW line of sight) and a 
density nH2=2000 cm~^ is ~3 pc. Gonsidering VLBI mea- 
surements of the size of t he continuum i mages (~2pc in 
the image plane at z=0.89. [jin et al1l2003f) . the high (near 
100%) covering factor of the absorbing clouds toward the 
SW line of sight, and the time variation of the absorption 
profiles, we conclude that the extent of the clouds should 
be of the order of a few parsecs. Hence, the bulk of the ab- 
sorption should be dominated by the cool component with 
21cin'^80 K. We thus estimate that the contribution from 
the warm component does not affect significantly our mea- 
surement of Tcmb , and our hypothesis of an homogeneous 
medium is reasonable. 

Although the derived value of the background radiation 
temperature is robust in our analysis, some neutral species 
like SiO and HNGO show slightly elevated rotation temper- 
atures, suggesting that they may require higher densities 
than the ions. This could be explained most easily if there 
is chemical stratification in the absorbing column. On the 
other hand, some of the crucial collision rates are poorly 
known, especially when the ortho form of H2 is the main 
collision partner. 

We also assume that the continuum source is fully cov- 
ered (/c=l in EqU]). Nevertheless, this is not a severe as- 
sumption for optically thin lines, and recent ALMA obser- 
vations of the ground-state line of ortho water, resolving 
both NE and SW images of the quasar, show that the SW 
continuum source covering factor is indeed close to unity 
(Muller et al., in prep.), at least at the time of the obser- 
vations in 2012. 

Finally, we have consider ed Hg molec ules only as col- 
lision partners. As shown by iLisztl (|2012[ ) for the Galactic 
interstellar medium, molecular ions with large dipole mo- 
ments like HGO"*" are most sensitive to coUisional excita- 
tion by electrons. Electron- impact begins to compete with 
neutral collisions when the fractional ionization 7i(e)/nH <^ 
10^^. The observed rotational temperatures of H^'^CO^ 
and HOC+ are both close to the expected Tcmb=5.14 K, 
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suggesting that the contribution from collisions with elec- 
trons is not significant. 

The robustness of the Tcmb measurement toward 
PKS 1830—211 could be improved by investigating the dif- 
ferent lines of sight to the lensed images of the background 
quasar. New observations when the NE absorption be- 
comes strong er again {e.g. as in the period 1999-2003, see 
iMuller fc Guelin ,2008h . with angular resolution sufficient 
to resolve the NE and SW images (separated by 1"), would 
provide an interesting independent Tcmb measurement. 

Additional transitions/molecules could be observed 
(within a short time interval to minimize time variations) 
to place tighter constraints on the physical conditions. 
However, and very importantly, collisional rates (with H2) 
are still lacking for a large number of molecules. This work 
is a good illustration of the importance and needs in astro- 
physical studi es of molecular sp ectroscopy and collisional 
data (see e.g. Ivan der Takll201ll for a short review on the 
field). 

We have shown that using molecular transitions as diag- 
nostics of physical conditions toward molecular absorbers 
is a powerful method to measure Tcmb- Unfortunately, 
there is only a small number of z>0 millimetre-wave 
band molecula r absorber s known to date (se e a review by 
ICombesI [200I and also ICurran et all 1201 ll ). Due to the 
brightness of the background continuum, the large col- 
umn density of absorbing gas, the large number of detected 
molecules, and the multiple independent lines of sight, the 
absorber toward PKS 1830—211 is a prime target for a pre- 
cise measurement of Tcmb- Another relatively well-suited 
target could be the z=0.68 absorber located toward the 
quasar B 0218-1-357 (|Wiklind fc Combes! 1199,51) . The phys- 
ical conditions have already been explored using transi- 
tions of ammonia jH cnkel et al. 2005) and formaldehyde 
(jJethava et aIll2007nZeiger fc'Darlindl2010[) . but a thor- 
ough investigation combining several molecules remains to 
be done for this source, in order to derive a precise mea- 
surement of Tcmb at z=0.68. A potentially severe difficulty 
is the small angular separation of ^0.3" between the two 
gravitationally lensed images of the quasar. 

In the (near) future, new z>0 molecular absorbers 
might be discovered (e.g. with the Atacama Large 
Millimeter-submillimeter Array or the Karl G. Jansky Very 
Large Array), providing more chances to test the Tcmb evo- 
lution with redshift using molecular excitation analysis. 

5. Summary and conclusions 

We have performed a multi-transition excitation analysis 
of various molecular absorption lines toward the radio-mm 
molecular absorber PKS 1830— 211, with the aim to pre- 
cisely determine the cosmic microwave background temper- 
ature at z=0.89. We use new 7 and 3 mm ATCA observa- 
tions obtained in 2011, co mplemented with 20 09 data from 
our ATCA 7 mm survey (jMuller et al.ll201lD . Our results 
can be summarized as follows: 

— The rotation temperatures of the observed molecules are 
ah close to ~5 K, comparable to the value Tcmb=5.14 K 
predicted from the adiabatic expansion of the Universe. 
This suggests that the molecular excitation is mostly 
radiatively coupled with the CMB, and that an accurate 
value of Tcmb can be derived. 



— We combine a Monte-Carlo Markov Chain approach 
with non-LTE radiative transfer predictions from the 
RADEX code to solve for the physical conditions of the 
absorbing gas, taking into account time variations of 
the quasar {i.e. possible changes of the sightline) be- 
tween the 2009 and 2011 observations. For each epoch, 
we assume that the absorbing screen is homogeneous, 
and characterized by a single kinetic temperature and 
H2 density. Under this (reasonable) assumption, and 
with the molecular data available, we determine a CMB 
temperature of 5.08±0.10 K (68% confidence level) at 
z=0.89, which is the most precise measurement of Tcmb 
at z>0 to date. 

— Combining this result with other measurements pub- 
lished in the literature, we refine the evolution of 
the CMB temperature as a function of redshift, 
T=Tox(l-Hz)(i-"), finding a=0.009±0.019, ie. a quasi- 
linear dependence to the redshift. As far as we can 
tell, this is consistent with adiabatic expansion of the 
Universe. 

— The absorbing gas has properties comparable to 
Galactic diffuse/translucent clouds, with a kinetic tem- 
perature of ~80 K and a H2 density of the order of 
10^ cm-3. 

— We report the detection of several new species toward 
PKS 1830-211: NH2CHO, ^"SiO, HCS+ and tentatively 
HOCO+. The weak lines of these species could be de- 
tected in 2011 data toward the SW component owing to 
the general increase in absorption depth of all species by 
a factor of two compared to our previous observations 
in 2009. 

This work emphasizes the extreme importance of accu- 
rate molecular collisional data for astrophysical studies. 
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Fig. 1. Spectra of the strong absorption of HCO+, HCN and HNC transitions observed in 2011 with ATCA (full black 
line). The dotted spectra correspond to PdBI observations with the narrow-band corr elator, i.e., a t better velocity 
resolution (see Table [T]). The J=l-0 spectra observed in our previous ATCA 7 mm survey (iMuUer et al.llMlD are shown 
in light grey. The intensity is normalized to the total continuum flux (NE plus SW images). 
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Fig. 2. Evolution of the column density (A^lte) as a function of the rotation temperature (Trot) for the transitions 
of different molecular species observed in 2011 f Table ETTj) . according to EqH) Each curve corresponds to a different 
transition and the dashed curves delimit the Icr uncertainties. The curve intersections give the values of A'coi and Trot- 
The vertical dotted line indicates Tcmb =5. 14 K. 
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Fig. 3. RADEX predictions of the rotation temperatures as a function of H2 density for the transitions of H^'^CO+, 
H^^CN, HC^^N, H NC, SiO, and H NCO hsted in Table EH and transitions of C-C3H2, I-C3H2, SO, and HC3N hsted 
in Table A.l bv Muller et al.ll201lL We assume Tcmb=5.14 K, Tkin=50 K (blue), and Tkin=100 K (red). The rotation 
temperature derived from the observations of two (or more) transitions is the average of the corresponding rotation 
temperatures curves. 
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Fig. 4. One and two dimensional lilceliliood distributions for the Metropolis-Hasting estimates of the non-LTE condition 
toward the SW component of the 2009 and 2011 campaign. Only a subset of the parameters is represented. The blue 
dashed lines correspond to the common range of temperatures for which collisional rates are available in the LAMDA 
database for the considered set of species. The normalized ID likelihood distributions {last row) are fitted by a Gaussian 
distribution (red continuous line). 
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Fig. A.l. Spectra from ATCA and PdBI observations. Tlie frequency scale is in the observer frame. The spectra 
normalized to the total flux density corresponding to the sum of the NE and SW images. 
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Fig. A. 2. Same as Fig lA.li but zoomed in for weak lines and with 2009/2010 averaged 7 mm spectrum overlaid (magenta). 
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Table A.l. Parameters and results of the global fit of the spectra for 
non-saturated lines detected toward the SW component in 2011 ATCA 
observations. 



Rest Freq. 


Transition 




Ei/kB 


Vo ' 


AV 


/ TdV " 
(10^3 km s-i) 


(MHz) 






(K) 


(km s ^) 


(km s^"'") 


87316.898 


C2H N=l-0 


1.00 


0.0 


-0.5 (0.1) 

_ t 


18.9 (0.1) 


28234 (280) 


174663.199 


C2H N=2-l 


2.00 


4.2 


_ t 


37574 (486) 


86339.922 


H"CN J=l-0 


1.00 


0.0 


-2.7 (0.1) 

„ t 


16.8 (0.2) 


4966 (69) 


172677.851 


Ri^CN J=2-l 


2.00 


4.1 


_ t 


6393(141) 


86054.966 


HC^^N J=l-0 


1.00 


0.0 


_ t 


_ t 


945(42) 


172107.957 


HC^^N J=2-l 


2.00 


4.1 


_ t 


_ t 


943(99) 


90663.568 


HNC J=l-0 


1.00 


0.0 


0.17 (0.03) 
_ t 


19.8 (0.1) 


43339 (1130) 


181324.758 


HNC J=2-l 


2.00 


4.4 


_ t 


45856 (1347) 


87090.825 


HN"C J=l-0 


1.00 


0.0 


_ t 


_ t 


1515 (50) 


174179.411 


HNi^C J=2-l 


2.00 


4.2 


_ t 


_ t 


1778 (86) 


88865.715 


Hi^NC J=l-0 


1.00 


0.0 


_ t 


_ t 


300 (50) 


86754.288 


Hi'^CO+ J=l-0 


1.00 


0.0 


-2.7 (0.1) 


15.6 (0.1) 


6089 (79) 


173506.700 


H"CO+ J=2-l 


2.00 


4.2 


_ t 


_ t 


8072 (153) 


85162 223 


HC^*0+ T— 1-0 


1.00 


0.0 


_ t 


_ t 


2581 (46) 


87057 535 


TTCi''0+ T— 1-0 

XiVJ V y — 1. yJ 


1.00 


0.0 


_ t 


_ t 


262 (40) 


1 741 1 1 fiQ 


Hri'^n+ T— 9 1 


9 no 


4 9 


_ t 


_ t 


9m /7n^ 




1 1 V 7 V. J — 1-U 


1 no 

1 .uu 


n n 

u.u 


— U.Z \^U.tJ J 




1 n^7 ((^oS 


1 7SQ79 n'^1 


nnp+ T— 9 1 

1 ivy V. J- — Zf-x 


9 on 

z.uu 


4 


_ t 


_ t 


1 98"? /q9^ 


8fifi7n 7fin 
OUU ( U. ( uu 


npo N,^^— ir^-, n^r, 1—^/9 1 /9 
1 1 V v7 K — ^01~'-J00 " — O/z-i/z 


1 on 

1 .uu 


u.u 




1 S 8 ('9 4') 
10.0 ^z.^^ 


000 \\jZj j 


1 72266 8'i3 




7 40 


9 9 


—4 1 ro 6"! 


9 8 n 41 


81 9 ?61 1 


89956 068 


\^ii3iNii2 I^Jjy ±10 -LQl 


6.00 


2.1 


—4.8 (1.5) 


12 3 ^3 2^ 


223 (hVi 


85338 893 


1^ v_^3ii2 -^12 J^Ol 


4.50 


2.3 


— 1.9 (0.1) 


17 8 ^0 2^ 


6973 f85^ 


84727 696 




96 


12.1 


_ t 


_ t 


944 ?QSl 


91958 726 


v_^i \ " — ■-'4: -'-4 


3.60 


123.1 


— 3 9 fO 5^1 


14.6 (1.1) 




91971 130 




12.80 


73.1 


_ t 


_ t 




91979.994 




8.40 


37.4 


_ t 


_ t 




91985.314 


CH3CN =5l-4l 


9.60 


16.0 


_ t 


_ t 




91987.088 


CH3CN J^f =5o-4o 


10.00 


8.8 


_ t 


_ t 




87925.237 


HNCO 4o4-3o3 


4.00 


6.3 


-4.5 (0.4) 


14.2 (1.0) 


845 (51) 


175843.695 


HNCO 8o8-7o7 


8.00 


29.5 


_ t 


_ t 


257 (81) 


86846.960 


SiO J=2-l 


2.00 


2.1 


-3.0 (0.1) 


UA (0.2) 


3612 (60) 


173688.310 


SiO J=4-3 


4.00 


12.5 


_ t 


_ t 


1944 (96) 


85759.199 


29SiO J==2-l 


2.00 


2.1 


_ t 


_ t 


541 (37) 


84746.170 


30SiO J=2-l 


2.00 


2.0 


__ t 


_ t 


287 (35) 


85347.890 


HCS+ J=2-l 


2.00 


2.0 


-1.5 (0.0) ^ 


15.0 (1.4) 


542 (45) 


85531.512 


HOCO+ 4o4-3o3 


4.00 


6.2 


0.5 (2.3) 


17.8 (5.5) 


181 (48) 


84542.329 


NH2CHO 4o4-3o3 


12.00 


6.1 


-4.4 (0.8) 


15.0 (0.0)^ 

_ t 


507 (52) 


87848.873 


NH2CHO 4i3-3i2 


11.23 


9.3 


_ t 


290 (42) 


178605.403 


SO Jk=54-43 


4.91 


15.9 


-3.4 (0.7) 


13.7 (1.6) 


843 (85) 


90979.023 


HC3N J=10-9 


10.00 


19.6 


-2.5 (0.7) 


14.6 (1.6) 


680 (66) 



Notes. " Line strength. ^ Energy of the lower level of the transition. Line centroid. FWHM linewidth. Integrated opacity. ' 
The parameter was tied with the precedent entry. ^ The value was fixed. For some species, the fit results might be slightly different 
(although within uncertainties) to the values given in Table [3l due to the different fitting procedure with tied parameters. 
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